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1. INTRODUCTION 

With the start-up of the LHC approaching, it is certainly expedient to make an 
overview on what we currently know on quarkonium production, both in proton- 
proton and heavy-ion collisions and on what we can expect from analyses to be 
carried at the LHC. Heavy-quarkonium production mechanism has always been - 
and still is - a subject of debate (for reviews see [1, 2, 3, 4, 5]). Heavy quarkonia 
have been often suggested as ideal probes in studies and analyses of complex 
phenomena. However, reality was later found to be much less simple than initially 
thought. A well known example is the suggestion to measure the suppression of 
J/ip production in heavy-ion collision as a smoking-gun signature of the creation 
of the quark gluon plasma (QGP) [6]. However, cold nuclear matter effects, such 
as shadowing, energy loss, absorption, etc. , were shown to play an important role 
and had to be considered in the interpretation of the experimental measurements. 
Furthermore, effects of the successive dissociations of higher-excited states which 
can decay into J/ip (■?/'', Xc) had to be taken into account, as well as more specific 
issues related to the description of the plasma itself. 

In fact, even in a much "cleaner" environment, such as in pp collisions, under- 
standing quarkonium production has been a challenge since the first measurements 
by the CDF Collaboration of the direct production of J/ip and ip' at ^/s = 1.8 
TeV [7, 8]. It is fair to say that, at present, a consistent theoretical picture that 
predicts both cross sections and the polarisation measurements for charmonium at 
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the Tevatron [9], along with the cross section from PHENIX at RHIC [10] is not 
available. For instance, the long-standing prediction of Non Relativistic QCD [11] 
(NRQCD) on the transverse polarisation of ip^s at high transverse momentum is 
not supported by the data. The most natural interpretation of such flagrant failure 
of NRQCD is that the charmonium system is too light for relativistic effects to be 
neglected. 

Indications that this might be indeed the case come from the agreement between 
theory and the available experimental data for T production in pp (and inclusive 
decays). In this case, relativistic corrections are expected to be less important 
and the leading state in the Fock expansion, i.e. the heavy-quark pair in a colour 
singlet ^Si to be dominant. The latest NLO predictions [12] which include some of 
the important NNLO a| corrections, show a satisfactory agreement with the data 
coming from the Tevatron [13, 14]. Once again, much is expected from polarisation 
measurements at the Tevatron and the LHC to conflrm that at least bottomonium 
predictions are well under control. 

In pA collisions, LHC measurements will certainly be of greatest importance to 
pinpoint the size of shadowing effects in the small x region which can also be studied 
via electromagnetic (aka. "ultraperipheral" ) AA collisions. Furthermore, data will 
allow to understand the absorption mechanisms at high-energy and subsequently 
to gain insights of the different formation time of the various heavy quarkonia. 

Finally, nucleus-nucleus collisions at the LHC will be the long-awaited ideal 
laboratory for the study of the QGP. Unprecedented temperatures will be reached; 
in conjunction with the high-luminosity beams, bottomonia will be at last promoted 
as a practical probe for the QGP formation. With the good news from higher-order 
Q CD-correction studies centered on T, we are certainly at the beginning of very 
exciting discovery years. 

2. EXPERIMENTAL CAPABILITIES 

2.1. ALICE (by Andry M. Rakotozafindrabe) 

ALICE is the LHC dedicated heavy-ion experiment. Its main physics goal is 
to study the properties of the hot and dense deconflned hadronic matter which is 
expected to be created during the relativistic heavy ion collisions. ALICE'S primary 
interest into the production rate of the heavy-quarkonium states lies into the fact 
that it can be used as a sensitive probe to the formation of the quark-gluon plasma. 
At temperatures above the quarkonium binding energy, the latter is foreseen to 
melt through colour screening, inducing a suppression of the production rates. The 
dissociation temperature T^ pattern would be: 

T,[^'] ^ T,[xc] < T,[T(35)] < T,[J/^] ^ T,[T{2S)] < T,[T{1S))] (1) 

which shows how quarkonium suppression can be used to estimate the temper- 
ature of the created QGP. Lower energy accelerators/colliders (SPS and RHIC) 
have indeed observed J/ijj suppression. For the T family, this will be only feasible 
at the LHC where the bb production cross-section is quite sizeable, and where 
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Trf[T(15')] can be reached. Moreover, at nominal luminosity, central PbPb collisions 
at the LHC {^/s~^ = 5.5 TeV) are expected to produce about a hundred of cc 
pairs, which substantially increases the regeneration probability of secondary 
J/^'s. This emphasises the interest in measuring the T(25') production rate, since 
it is expected to benefit from a small regeneration probability (about only five 
bb pairs are expected in central PbPb collisions). An important remark is that 
in-medium effects on charmonia can be studied once the feed-down from B decays 
is properly subtracted {B -^ J jip + X, expected to account for about 20% of the 
total J/ip yield if no cold nuclear matter effect is considered). Obviously, results 
obtained in AA collisions must be benchmarked against the ones from pp col- 
lisions, but also pA collisions in order to get the baseline for the cold nuclear effects. 

The LHC program plans pp running at ^/s = 14 TeV for 8 months per year (10^ s 
effective time) followed by ion running for one month per year (10^ s). ALICE can 
participate to pp running, but at a maximum allowed luminosity of 5 10^° cm~^ s~^. 
The first heavy-ion run will be PbPb collisions at y/s^ = 5.5 TeV at a luminosity 
of 510^^ cm~^ s~^, corresponding to 1/20**^ of the design luminosity. One or two 
years of light-ion collisions {ArAr, with a luminosity up to 10^*^ cm~^ s"^)) and one 
year of pPb collisions at y/s^ = 8.8 TeV are also planned. 

ALICE [15, 16] can detect quarkonia in the dielectron channel at central rapidity 
{\ri\ < 0.9), and in the dimuon channel at forward rapidity (—4 <r] < —2.5). In these 
rapidity ranges, the corresponding probed Bjorken-x approximately goes from 10~^ 
to values as low as 10~^, providing additional insight into the PDFs and their 
modification in the nuclear environment. The central barrel and the muon arm are 
equipped with dedicated triggers on the individual electron and muon transverse 
momentum. In order to select a lepton pair from quarkonium decays, a 1ow-Pt 
trigger cut is applied on individual leptons to remove most background sources: 
Pt >3 GeV/c for single electrons and Pt > 1 (2) GeV/c for single muons from 
charmonia (bottomonia). As a consequence, it prevents the detection of charmonia 
with Pt < 5 GeV/c in the dielectron channel, whereas the charmonia can be 
detected down to very low P^ (about a hundred of MeV/c) in the dimuon channel. 
The high Pt reach is expected to be 10 (20) GeV/c for the J/ip into dielectrons 
(dimuons) for a PbPb run of one month at nominal luminosity. In both channels, 
the expected mass resolution of about 90 MeV/c^ will be sufficient to resolve all 
the T states. The expected mass resolution for the J/ip is about 30 (70) MeV/c^ 
in the dielectron (dimuon) channel. The central barrel has excellent secondary 
vertexing capabilities, combined with particle identification. Therefore, prompt 
and secondary J/ip, from B decays, can be distinguished at central rapidities via 
a displaced vertex measurement. At forward rapidity, this technique is not usable: 
the prompt J/ip yield has to be determined indirectly, by subtracting from the 
measured yield the one expected from B decay. The latter is inferred from the 
single- muon P^-spectra measurement with the cut P^ ^1-5 GeV/c applied to all 
reconstructed muons to maximise the beauty-signal significance. A fit technique is 
then applied to extract the Pt distribution of the muons from B decays [16]. Last, 
but not least, ALICE will be able to measure the J/ip polarisation, both in pp and 
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PbPb collisions. 



2.2. ATLAS and CMS (by Aafke C. Kraan and David d'Enterria) 

The ATLAS [17] and CMS [18] experiments at the LHC are general purpose 
detectors designed to explore the physics at the TeV energy scale. The primary 
goals of the experiments are to reveal the electroweak symmetry breaking mech- 
anism and provide evidence of physics beyond the Standard Model (SM) in 
proton-proton collisions at ^/s = 14 TeV. The two experiments are obviously 
extremely well suited to carry out Quantum-Chromodynamics studies in both pp 
and PbPb collisions. 

The total cross section for J/ip and T production at pp collisions at 14 TeV are 
expected to be around 0.4 mb and 7 /ib, respectively, and for the higher mass states 
an order of magnitude lower. The J/ip and T measurements in ATLAS and CMS 
focuse normally in the dimuon decay channel {J/tp^T — > /i"'"/i~) with branching 
ratios of J /ip —>■ fi~^fJ,~ and T — > /i"*"//" are 5.98% and 2.48% respectively. On the 
one hand, the high centre-of-mass energies (up to 14 TeV) and luminosities (up 
to 10^^ cm^^s^^) as well as the large-acceptance muon spectrometers (|?7| <2.5, 
full azimuth) available in the ATLAS and CMS experiments will allow, unlike 
the ALICE and LHCb cases, for quarkonium measurements at the LHC up to 
very large transverse momenta. On the other hand, the strong magnetic field 
and/or the large material budget results in lower reconstruction efficiencies for 
small transverse momentum muons (usually below Pt ^ 3 GeV/c). Thus, the 
sensitivity for quarkonium measurements at low momentum is inferior than that 
of the Tevatron or ALICE and LHCb experiments. Only startup conditions, when 
the luminosities and muon trigger thresholds will be low, will allow for some lower 
transverse momentum measurements. 

Besides physics measurements, during detector commissioning and startup, 
dimuon invariant masses of quarkonium states will provide extremely useful can- 
dles for detector calibration and alignment. Several physics analyses are planned 
at both ATLAS and CMS including: 

• Differential inclusive cross section for J/ip, T, ip{2S), T(2S'), T(3S'), as well 
as the Xc, x\i xl states in pp, pA and A A colhsions. 

• Polarisation measurements of these states. 

• More exclusive measurements aimed at understanding the underlying QQ 
production mechanisms by looking, e.g., at the associated hadronic activity. 

In Table 1 some relevant parameters for J/ip^s and T's are given for the ATLAS 
and CMS experiments [17, 18, 19, 20]. For the heavy-ion running, the performances 
of ATLAS [21] and CMS [22] are very similar to the proton-proton ones. 
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TABLE 1. 

at 14 TcV. 



Basic J /ip and T reconstruction performances for ATLAS and CMS in pp collisions 



1 


ATLAS 




CMS 


II 


QQ trigger threshold 


P^^ > 6 GeV/c, P^' 


> 4 GeV/c, 


P^^ > 3 GeV/c, P^^ 


> 3 GeV/c 


J/tp mass resolution 


-50 MeV/c2 




-30 MeV/c2 




iV-™ , in 10 pb-i 


2x10^ 




3x10^ 




T mass resolution 


~170 MeV/c2 




~95 MeV/c2 




AT-- ^ i^ 10 pb-i 


5x10^ 




1x10^ 





2.3. LHCb (by David d'Enterria) 

The LHCb experiment [23] at the LHC is mainly focused on the search of 
possible signals of new physics in CP-violation and rare decays processes in the 
heavy-quark sector of the Standard Model (SM). For this purpose, the experiment 
has been equipped with arguably the best capabilities for the detection of h and 
c quarks produced at forward rapidities in proton-proton collisions at the LHC. 
Indeed, the LHCb detector - with a single-arm configuration - has excellent and 
varied particle detection and identification capabilities in the forward hemisphere. 
In particular, muons and electrons can be well measured in the pseudo-rapidity 
range 1.8 < t] < 4.9. For comparison, in the same rj range ATLAS and CMS 
can only reconstruct electrons (in a reduced 2 < r^ < 3 range), whereas ALICE 
can only measure muons (in a slightly more reduced acceptance: 2.5 < r/ < 4, 
and only for proton-proton luminosities, 10^*^ cm~^s~^, 100 times lower than those 
available for LHCb). Those characteristics make of LHCb an excellent apparatus 
to measure forward quarkonium production cross-sections and polarisation via 
J/i'yT -^ e+e~, /i"*"//", including the excited states {ip{2S), T(2S), T(3S')). 



3. OPEN ISSUES 

3.1. pp collisions (by Andry M. Rakotozafindrabe and J. P. Lansberg) 

The underlying theory for direct and prompt ip is still under intense debate [1, 2]. 
Via the colour octet (CO) mechanism, NRQCD factorisation [11] has been success- 
ful to explain some features of the charmonium hadroproduction. As illustrated 
by the comparisons to CDF measurements in pp [7, 8] or to PHENIX old mea- 
surements in pp [24, 25]), for -Pr ^ 5 GeV/c, it provides a good description of 
the Pt- differential cross-section for the direct J/tp and tp' , the cross-section being 
dominated by the gluon fragmentation into a colour-octet ^5*1 state. The latter 
mechanism leads to transversally polarised J/ip and ip'. 

However, this is not seen by the CDF experiment [9] which measured a slight 
longitudinal polarisation for both the prompt J/ip and direct ip' yield. It is worth 
noting here that the feed-down from Xc can influence signiflcantly the polarisation 
of the prompt J/ip yield - this was taken into account in the NRQCD-based predic- 
tions [26]. Moreover, the recent preliminary result from PHENIX [27] indicates a 
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polarisation compatible with zero for the total J /ip production at forward rapidity 
(1.2 < \y\ < 2.2), but with large uncertainties. 

It is therefore not surprising to observe a renewed interest in improving the 
present predictions for the colour singlet (CS) contribution, by computing the 
higher-order QCD (see section 4.1) corrections in as, or by "softening" some 
of the basic assumptions of the common approaches, as done in [28] via the 
consideration of the s-channel cut contribution. On the one hand, NLO [29, 30] and 
part of NNLO corrections [12] significantly enhance the quarkonium yields^. In the 
T(15') and T{3S) cases, the dominant NNLO corrections to the colour singlet [12] 
suffice to successfully describe the measured Pr-differential cross-section of the 
direct yield [13, 14]. The polarisation predictions for the latter cases seem quite 
encouraging considering CDF [13] and D0 [32] measurements. Those corrections 
could be though still unable to bring agreement with the measured P^- differential 
cross-section of the direct J/ip, but dedicated further studies are needed. 

On the other hand, by including s-channel cut contributions [28] to the usual 
production CS production mechanism one can reproduce the Pr-spectra up to 
intermediate values of the J/ip^s Pt, both at the Tevatron and at RHIC, and provide 
mostly longitudinally polarised J/ip at the Tevatron. However, as expected [28], 
this approach underestimates the cross-section at large values of Pt- 

In summary, the theoretical status sounds clearer for the bottomonia than for 
the charmonium production processes. Additional tests are undoubtedly needed 
beyond the mere measurements of inclusive cross section and polarisation at 
the LHC. For instance, the hadroproduction of J/ip or T with a heavy-quark 
pair [30, 33] appears to be a new valuable tool to separately probe the CS 
contribution, at least dominant at low-Py (below 15 GeV/c), as well as the study 
of the hadronic activity around the quarkonium (see section 5). 



3.2. pA collisions (by Andry M. Rakotozafindrabe) 

The interest of pA collisions is based on the possibility they open up to evaluate 
both the initial and final-state effects on heavy-quarkonium production in cold 
nuclear matter (CNM). Such baseline is mandatory to be able to draw conclusions 
about any further effects due to QGP formation in AA collisions. In the following, 
we will address a non-exhaustive list^ of these CNM effects. 

Initial-state effects. Heavy-quarkonium production mainly proceeds through 
gluon fusion at relativistic ion colliders. Therefore, the nuclear shadowing of initial 



^ This sounds like a confirmation of the study [31] which dealt in a simplified way with NNLO 
corrections assimilable to LO BFKL contributions. However, this study could not provide a 
prediction for the Pt dependence: it only predicted an enhancement of some NNLO corrections 
for large s. 

For instance, let us mention a potential further significant effect due to charm-quark shadowing 
as computed in [34] in a light-cone Green function formalism. 
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gluons has been extensively investigated (see [35] for a recent review), together 
with its consequences on the charmonium production. 

On the experimental side, the gluon nPDF"^ is loosely constrained at small values 
of Bjorken-x: 

• On the one hand, the processes used - Deep-Inelastic scattering (DIS) and 
Drell-Yan - are mostly sensitive to the quark and antiquark densities. There- 
fore, the gluon density is indirectly constrained, either via the nucleon- 
structure-function deviations from the Bjorken scaling, caused by gluon ra- 
diation, or via sum-rules (conservation of the nucleon momentum distributed 
among all partons). 

• On the other hand, there is no nuclear DIS data below x < 5 10~^ at pertur- 
bative values of the momentum transfer Q^ > Aqqj-,, required for the validity 
of the DGLAP equations used to predict the evolution of nPDF with Q^ . 

As a result, the extracted parametrisations of the ratios nPDF/PDF have large 
uncertainties at low-x: typically, for the gluon shadowing in Pb, the LO parametri- 
sations EKS98 [36] and EPS08 [37] give values of nPDF/PDF that differ by about 
a factor of ten at a; = 10~^ and Q^ = 1.69 GeV^. EPS08 notably includes addi- 
tional constraints from high-P^ {Pt > 2 GeV/c) hadron production measured by 
the BRAHMS experiment [38] at forward rapidities in dAu collisions at RHIC top 
energy. These data mainly probe gluons with a; > 5 10~^ in the gold nucleus and 
suggest a much stronger shadowing. The discrepancy is even larger when these LO 
parametrisations are compared to the NLO ones, either nDSG [39] or HKN07 [40]. 
All these uncertainties preclude yet reliable predictions of heavy-quarkonium pro- 
duction at the LHC, dominated by low-x gluons (see Section 2.1). 

A workaround to the use of these parametrisations can be found in approaches 
that try to describe the shadowing in a formal way: 

• The underlying physical mechanism is thought to be multiple-scattering (or 
multi-pomeron exchanges) with initial interactions between the pomerons, and 
the calculations are made in the Glauber-Gribov framework [41]. However, 
these models usually have a narrower validity range (limited to the low-x 
region [42]) since they were designed to describe the coherence effects that 
lead to the depletion of the nuclear structure function. Indeed, at intermediate 
values of x, the amount of anti-shadowing appears to be smaller in these 
approaches than in the aforementioned parametrisations (for instance, see the 
comparisons reported in [43] for the J/ip shadowing in dAu at RHIC). 

• The Colour Glass Condensate (CGC) is an effective theory which describes the 
behaviour of the small-x components of the hadronic wavefunction in QCD 
(see [44] for a recent review). Hence, it can be used to study the high energy 
scattering in QCD, namely the initial stages of heavy-ion collisions. The theory 
is characterised by a saturation scale Q^: at any Q^ below this scale, the rapid 
rise of the gluon density at small- a; slows down to a logarithmic rate, due to 



nPDF stands for the parton density within a bound nucleon. 
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a growing number of gluon-gluon fusions. It is of no doubt that pA collisions 
at the LHC will provide crucial tests of the CGC framework, since they will 
allow to deeply probe the saturation region. At RHIC, the understanding of the 
initial effects on charmonium production in the CGC framework is still a work 
in progress. A first step was the calculation [45] of the open charm production, 
which is found suppressed at forward rapidity at RHIC. Predictions for the 
J/ip are under way: for rapidities y >0, a rather qualitative agreement with 
PHENIX dAu data is obtained in [46] ^ for the ?/-dependence and for the 
centrality dependence ^. 

To summarise, the amount of shadowing crucially depends on x. By taking the 
J/ip transverse momentum Pt into account when evaluating Xi, X2 (and Q"^), the 
influence of Pt on the shadowing is investigated in [43, 51] at RHIC. There is an 
on-going debate on the way the J/ip acquires its Pt: either (i) the initial gluons 
carry an intrinsic transverse momentum and the latter is subsequently transferred 
to the J/ip (2^1 process), or (ii) the Pt comes from the emission of a recoiling 
outgoing gluon (2 — > 2 process). For the latter process, the authors consider the 
partonic cross-section given in [28] which satisfactorily describes the J/ip Pt- 
spectrum down to Pt ^ at RHIC. On the average, initial gluons involved in 
these processes originate from different x-regions, hence resulting in quite different 
shadowing effects [51]. The present uncertainties on PHENIX dAu data [49] do 
not allow to discriminate these scenarios, but forthcoming improvements will be 
obtained from the recent data taking (with at least a factor of 30 in statistics). 

Additional initial-state effects are the initial-partonic multiple scattering and 
the related parton energy loss. It is believed that the observed broadening of (P^) 
- the mean value of the J/^'s transverse momentum squared - from pp collisions 
to pA collisions (with increasing A) is due to such multiple-scattering (the so- 
called "Cronin effect"). This effect is usually described as a random walk of the 
initial-projectile parton within the target nucleus (see e.g. [52, 53]): the resulting 
(Py) proportionally increases with the amount of scattering centers, characterised 
by the length L of nuclear matter traversed. Simple linear fits to (P^) vs L can 
indeed account for the pA measurements done at SPS. All pA and AA results at 
SPS, but the preliminary result reported by NA60 in pA at 158 AGeV, exhibit 
the same slope [54]. At RHIC, the (P^) measured in dAu [49] suffers from large 
uncertainties, but is compatible with a moderate broadening. A linear fit done to 
all data points available at RHIC {pp, dAu, CuCu and AuAu) is reported in [55]: 
the slope is compatible with zero at mid rapidity and with some broadening at 
forward rapidity. Interestingly, within the large uncertainties, the slope seen at 
forward rapidity at RHIC is compatible with the one at SPS (the comparison can 
be made between the slopes quoted in [55] and in the slides [56]). 



^ Note that the data to theory comparison made in this article should be updated with the newly 
published PHENIX dAu results from the re-analysis described in [49] . 

This work is being extended in [50] in order to describe the y-dependence of the peripheral 
AuAu collisions at RHIC. The first results seem quite promising. 



Perspectives on heavy-quarkonium production at the LHC J. P. Lansberg et al. 



Final-state effects. The so-called "nuclear absorption" has been extensively 
investigated for the charmonia. It reflects the break-up of correlated cc pairs due to 
inelastic scattering with the remaining nucleons from the incident cold nuclei. As 
we shall see below, the underlying mechanism is still unclear. Moreover, its physical 
picture may change with energy, as pointed in [57] and the transition would occur 
at RHIC energies: 

• At low energy, there is a longitudinally-ordered scattering of the heavy-quark 
pair. It results in an attenuation factor exp[—aahsPQL{b)], where cxabs is the 
effective break-up cross-section, po = 0.17 nucleonfm"^ is the nuclear density 
and L{h) is the length traversed by the heavy-quark pair in the nuclear 
matter at a given impact parameter h. The break-up cross-section is not 
calculable from flrst-principle QCD^. Presently, it is a free parameter in the 
models, such as in [49] where its value is obtained by fltting the data with 
a given nuclear shadowing model and an unknown additional absorption. 
Many considerations are hiding behind the effective value of Uabs- It has 
been argued [58, 59, 60, 61] that the value of the break-up cross-section 
should depend both on the collision energy and on the colour state of the 
created cc. With increasing collision energy, the pair will hadronise from 
inside to outside the nucleus. Since a hadronised cc is much more robust 
than the pre-resonance, the effective value of Uabs will decrease with energy. A 
colour singlet pair has a smaller size and hence a shorter hadronisation time 
than a colour octet pair, so the corresponding break-up cross-section should 
be smaller. Moreover, if the direct J/ij) and ip' are believed to be created 
perturbatively in the same cc octet state, then they should suffer the same 
amount of nuclear absorption, unless the hadronised states are also broken-up 
by the scattering off nucleons. A quite rich pA program was developed at SPS, 
where the CNM effects can be described with the nuclear absorption only: the 
latest values from the NA50 experiment are [62] (Tahs{J/'4') = 4.2 ±0.5 mb and 
c"abs('0') = 7.7±0.9 mb. However, the quoted value of crabs(^/V') is for the total 
J/tp production: it accounts for the absorption of the direct J/ip as well as 
of the higher mass cc states that would have decayed in J/ip otherwise. It 
is also worth remembering that the shadowing effect was "forgotten" when 
evaluating the quoted values of (Tabs- When properly taking into account 
the anti-shadowing at SPS as in [61], the estimated break-up cross-sections 
are signiflcantly larger (around 7 mb for the J/ip and 11 mb for the ip'). 
At RHIC, the published values of the effective break-up cross-section in 
dAu colhsions [49] are crabs(^/^) = 2. 8+1:4(2. 2+};5) mb for the two different 
nPDF/PDF parametrisations used to evaluate the shadowing. However, the 
shadowing used in [49] is obtained without taking into account the J/ip^s Pt- 
When doing so, in the 2 — >■ 2 process, the resulting shadowing is quite different, 
leading to a different value of the extracted break-up cross-section [51], which 



^ It is nevertheless possible to extract its value from photoproduction data 7A^ -^ ipN (see 
e.g. [47, 48]). 

Perspectives on heavy-quarkoniuni production at the LHC J. P. Lansberg et al. 10 



happens to be closer to the pubhshed NA50 value. 

At high energy, the above space-time picture of the heavy-quarkonium pro- 
duction in pA should not hold any more. As argued in [57, 63], the heavy 
state in the projectile will rather undergo a coherent scattering off the target 
nucleons. The conventional treatment of nuclear absorption is not valid any 
more and the rigorous method reported in [64] has to be used. As a result, (Tabs 
asymptotically tends to zero at high energies. The required threshold in ^/s 



NN 



for this coherent scattering picture to be relevant is a function of x+, the lon- 
gitudinal momentum fraction of the heavy system. The transition happens at 
RHIC energy for the J/ip produced at mid rapidity. Consequently, the break- 
up cross-section is rapidity-dependent at RHIC. In [57, 63], the authors use 
the method [64] within a Glauber-Gribov description of shadowing to predict 
the CNM effects on the J/ip production at various y/s~^, from 39 GeV to the 
LHC energies. They report a fair agreement with E866, E772 and PHENIX 
p{d)A data, both for the rapidity-dependence and the x^-dependence. 

The J/ip puzzle. As can be seen from the above discussion, disentangling the 
various CNM effects at work for the J/tp production in pA is not an easy task. So, 
in order to set a proper baseline for the QGP search, how can we safely extrapolate 
the CNM effects from pA to AA ? The current solution adopted by PHENIX [49] 
is to use the data-driven method inspired by [65] to predict the expected CNM 
effects in AuAu from the available dAu data. This approach has the advantage of 
a proper propagation of the experimental errors. The quite annoying drawback is 
that it is unable to derive any specie- nor energy- dependent extrapolation of the 
CNM effects. 

To improve this situation, we clearly need more and much precise data. At RHIC, 
the analysis of the latest high statistics dAu run is well under way. At the LHC, a 
rich pA program will be crucial. 

In the meantime, a key test to the models would be to reproduce the observed 
scaling of the a exponent in o"p(d)A = CTpp^", with xp for the J/ip production in 
pA at various ^/s^. At least three different approaches tried to cope with this 
observation. A first tentative is done in [59, 60], with a combination of five different 
types of CNM effects. More recently, with the help of the CGC approach and a 
simple (single- valued) (Tabs, the authors of [46] obtained a result with a reasonable 
agreement to the observed scaling. In the previous paragraph, we already mentioned 
the positive result achieved in [57] by the use of a coherent scattering picture of the 
nuclear absorption and a Gribov-Glauber framework to describe the shadowing. 



3.3. AA collisions (by Joseph Cugnon) 

The interest of heavy-quarkonium production in heavy-ion collisions is twofold. 
First, these collisions may offer different mechanisms of quarkonium production 
compared to pp collisions. Second, and more importantly, this production can 
be used as a sensitive probe to the supposedly formed quark-gluon plasma in 
the course of these collisions. There is a general consensus that, in the SPS 
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and RHIC conditions, the quarkonium states (mainly J /ip) are formed during 
the first nucleon-nucleon (NN) colhsions as in free space and that they are 
immersed afterwards in the quark-gluon plasma, where they can be dissolved due 
the screening of the quark-antiquark interaction in the plasma. In these conditions 
the plasma is rather cold, few charmed quarks are thermally produced and the 
final J/tp yield results from the absorption of the early produced charmonium 
states by the plasma. The situation will be somehow different at the LHC The 
temperature of the plasma will be substantially larger, free charmed quarks will 
be more numerous. Furthermore, charm production will result also from the decay 
of bottom, which will be produced sizeably at y/s^ = 5.5 TeV. In addition to 
the usual formation in first NN collisions, charmonium states will be produced 
also through the decay of i?-mesons and through the fusion of cc pairs, either 
thermally produced directly or through the decays of free 6-quarks (the so-called 
regeneration process). In relation to the use of quarkonia as probes of the plasma, 
these additional production processes should be clarified both experimentally and 
theoretically. First the side-feeding in the propagation of excited states in the 
plasma should be evaluated correctly. It should be reminded that the side-feeding of 
J/^|J is not yet satisfactorily clarified in the previous experiments [66]. As indicated 
in Section 2.1, this problem may be circumvented by using the T family as probes 
of the plasma. Nevertheless, the reliable evaluation of the feeding of J/ip by higher 
mass resonances in the plasma seems to be a necessary step on the analysis of 
future experiments. Furthermore, if one is interested in the energy or Pt spectrum 
of produced quarkonia, the study of the propagation of heavy quarks is highly 
desired. Progress have been made recently [67], with the use of Fokker-Planck 
equations (in relation to jet quenching), but the application of this approach is 
just beginning [68]. 

The discussion of our present understanding of quarkonium production in terms 
of known properties of these objects and the supposed properties of the plasma 
is postponed to Section 6. It probably requires a reliable description of the space- 
time evolution of the plasma and a good understanding of the interaction between 
a quarkonium and the surrounding medium. In the past, this problem has been 
tentatively circumvented by focusing on variables that are hopefully not sensitive 
to the detail of this evolution. The most popular example is provided by the famous 
plots of Raa versus L, the estimated average path of the quarkonium inside the 
plasma and admittedly a very crude variable. This has obvious limitations and 
does not presently provide a coherent view of the existing experimental data [69] . 
Therefore, an important effort in order to improve the theoretical description of 
the evolving plasma is certainly required. 

Concerning the coupling of the quarkonia to the plasma, new lattice-QCD 
developments have taken place recently. There are more and more indications that 
quarkonia may survive in the plasma up to temperatures of the order of two times 
the critical temperature Tc [70] . Whether the origin of this property results from 
unsuspected aspects of the screening properties of the plasma [71] or from other 
sources is beyond the scope of this overview. As a consequence, the situation at the 
LHC will be more comfortable than before, since the temperature of the plasma 
will be higher. In any case, a good theoretical understanding of the coupling is 
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necessary. It is not granted that this couphng is describable in terms of cross 
sections for quarkonium-gluon colhsions. If it is the case, recent progresses have 
been made in this field [72] , ahhough the variation of the cross section with energy 
seems to pose problems [66]. 

Finally, the diagnostic of the plasma will be possible only if the other probes 
are also well understood: high-Pi- photons, dileptons (coming from low mass reso- 
nances), etc. The first one is sometimes considered as the ideal probe. However, the 
present situation is not clear. There are inconsistencies between AuAu and CuCu 
data [66]. The second one seems to be more distorted than mesonic resonances by 
medium effects. Further theoretical work is necessary. 

4. RECENT THEORETICAL ADVANCES 

4.1. QCD corrections (by Jean-Philippe Lansberg) 

Recently, substantial progress has been achieved in the computation of higher- 
order QCD corrections to the hard amplitudes of quarkonium-production processes. 
The first NLO calculation to date was centered on unpolarised photoproduction of 
if) [73] via a colour-singlet (CS) state [74] (that is, LO in v for NRQCD [11]) more 
than ten years ago now. Later on, NLO corrections were computed for direct 77 
collisions [75, 76] where it had been shown [77] previously that the CS contribution 
alone was not able to correctly reproduce the measured rates by DELPHI [78]. 

At the LHC and the Tevatron, ip and T production proceeds most uniquely via 
gluon-fusion process. The corresponding cross section at NLO (q;| for hadropro- 
duction processes) are significantly more complicated to compute and became only 
available one year ago [29, 30]. Those results were recently confirmed in [79, 80]. 
In the latter papers, the polarisation information was kept and the observable a 
was also computed. It is important to stress that for ifj and T production the CS 
yields predicted at the NLO accuracy are still clearly below the experimental data. 
In this respect, the predictions for the polarisation at this order cannot be usefully 
compared to the data. 

Aside from hadroproduction regime, NLO corrections have also recently been 
computed for two J/ ^-production observables at the 5-factories: J /ip + cc [81] 
and J/i) + y]c [82]. 

The common feature of these calculations is the significant size of the NLO 
corrections, in particular for large transverse momenta Pt of the quarkonia when 
this information is available. In 7p an pp collisions, QCD corrections to the CS 
production open new channels with a different behaviour in Pt which indeed raise 
substantially the cross section in the large- Pt region. In general, the CS prediction 
can thus be brought considerably closer to the data, although agreement is only 
reached at NLO in the photoproduction case [73]. As of today, only the full colour- 
octet (CO) contributions to direct 77 collisions have been evaluated at NLO for 

Pt > [75, 76]. Very recently, CO contributions from S waves {^Sq and '^S{ ) have 
become available [83] for hadroproduction, but their impact on phenomenology has 
not been fully studied yet. Since NLO corrections do not affect significantly the 
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Pt dependence, a quick assessment can be obtained via the K factors, the ratios 
of NLO to LO predictions. The K factors of the cross section at the Tevatron 
are about 1.2 for the ^Sq state and 1.1 for the ^S*} one (at the LHC, they are 
both about 0.8). This entails that the value of the CO Long Distance Matrix 
Elements (LDMEs) fit to the Tevatron data at LO {0{^Sf^)) ~ 0.0012 GeV^ and 

(0(^4^')) ^ 0.0045 GeV^ [84] would be at most reduced by 15%. In this respect, 
the NLO corrections to the octets do not improve the universality of the matrix 
elements when the idea of the dominance of the CO transitions is confronted to 
the data on photoproduction from HERA. 

In [83], the authors made a fit of the CO LDMEs on prompt data which is 
not easily comparable to the previous works since it includes feeddown from ip' 
and Xc- It is however worthwhile to note that they had to abandon [in the fit] 
the experimental data with Pt < 6 GeV/c, since it is not possible to obtain a 
satisfactory Pt distribution in terms of a unique (O^) value. This emphasises the 
need for more work dedicated to the description of the low-P^ region. Last but not 
least, the polarisation from CO transitions appears not to be modified at NLO with 
respect to LO result, thus confirming the flagrant discrepancy between the NRQCD 
predictions for the polarisation of the J/ip and the experimental measurements 
from the CDF collaboration [9]. 

On the bottomonium side, the situation seems less problematic. We have now at 
our disposal the CS cross section at NLO including a dominant subset of NNLO 
corrections at q;| (namely the associated production with 3 light partons) for inclu- 
sive T hadropro duct ion [12] at mid and large Pt- The rate obtained by including 
only the CS channels are in substantial agreement with the experimental measure- 
ments of the cross section from the Tevatron [13, 14] . Concerning the polarisation, 
the direct yield is predicted to be mostly longitudinal. The experimental data be- 
ing centered on prompt yield [13, 32], we would need first to gain some insights 
on NLO corrections to P-wave production at Pt > to draw further conclusions. 
Yet, since the yield from P-wave feeddown is likely to give transversely polarised 
T, the trend is more than encouraging. 

4.2. Automated generation of quarkonium amplitudes in 

JNRQCD (by Pierre Artoisenet) 

The computation of heavy-quarkonium cross section within Non-Relativistic 
QCD [1 1] takes advantage of the small relative velocity v inside the quarkonium 
state to factorise in a consistent way perturbative high-energy effects (linked to the 
heavy-quark-pair production) from non-perturbative low-energy effects (linked to 
the evolution of the heavy-quark pair into a quarkonium state). This factorisation 
is performed order by order in as (the strong coupling constant) and in v, and is 
controlled by a factorisation scale A. As a result, differential cross sections read 

daiQ) = J^d^A (QQin)) {O^in))^ (2) 
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where n specifies the quantum numbers of the intermediate heavy-quark pair QQ. 
The factors da^ {QQin)) are called the short- distance coefficients, and can be 
computed perturbatively in as for a given process. The Long-Distance Matrix 
Elements (LDMEs) {0^{n))tv encode the soft evolution of the heavy-quark pair 
into a quarkonium state. They are universal, i.e., they do not depend on the details 
of the creation of the heavy-quark pair. 

In order to predict the cross section for a given process with a quarkonium 
in the final state, one has to compute the short distance coefficients at a given 
accuracy in as and in v (which limits the number of transitions n in the sum in 
Eq. (2)). Even for the tree-level amplitudes, a calculation by hand may be lengthy 
and in general error-prone, depending on the parton multiplicity in the final state. 
For the purpose of a fast, easy-to-handle and reliable computation of tree- level 
quarkonium amplitudes, a new MadGraph-based implementation [85] (MadOnia) 
has been developed recently. The user may require an S- or a P-wave, colour- 
singlet (CS) or colour octet (CO) intermediate state, for any process attainable in 
MadGraph for the open-quark production. The code then generates automatically 
the related squared amplitude at leading order in f , which can then be interfaced 
with a phase-space generator to obtain the short- distance cross section. For 5*- 
wave state production, the algorithm can also be extended to compute relativistic 
corrections to the Born- level cross section. For ^^i states, the decay into leptons 
can be included, thus giving an easy handle onto polarisation studies. Another 
capability of the code is the generation of amplitudes involving heavy quarkonium 
with mixed fiavors, such as the B^. 

Several applications of the code have already been reported [12, 30, 33, 85]. 
One example is the analysis of the associated production of a J /ip or a T plus a 
heavy-quark pair of the same fiavor at the Tevatron. The prediction of differential 
cross sections as well as polarisation observables, both for CS and CO transitions, is 
straightforward. One can directly show, for example, that the CS fragmentation ap- 
proximation, which was used to estimate the J/il) + cc production at the Tevatron, 
appears to underestimate the yield by a large factor in the region Pt < 20 GeV/c. 

Beside offering the possibility to check very efficiently the tree-level contribution 
of a large set of quarkonium production processes, the code is also embedded in a 
larger project aimed to serve the connection between theory and experiments. At 
the LHC, J/ip or T hadro-production followed by their leptonic decay will offer 
a very clean signature to calibrate the detectors as well as to probe new physics 
effects. Given the large number of such events, one can even hope to reconstruct 
more exclusive final states, such as T + 2 6-jets, and hence enlarge the number 
of measured observables to be compared with theoretical predictions. On the one 
hand, such measurements rely on theoretical assumptions to establish the criteria 
for the selection of the signal and to estimate the cut efficiency. On the other hand, 
theoretical predictions make use of the existing experimental data to quantify non- 
perturbative low-energy effects, including the values of the LDMEs in the case of 
quarkonium production. 

The flow of information between theory and experiment is made easier by the 
use of Monte Carlo tools. The new generation of these tools operate in two steps. 
First they produce parton-level events according to a hard scattering probability 
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which can be computed perturbatively in as- The events are then passed through 
a code that generates the parton shower and turns the partons into hadrons. 
Eventually one can use a detector simulator to smear the information according to 
the resolution of the detector, such that events are as close as possible to real data. 
For quarkonium production within NRQCD, the relative abundance of parton-level 
events is controlled by the short distance coefficients appearing in Eq. (2), which 
can be computed by MadOnia. The algorithm is being currently promoted to an 
event generator. Spin correlation and colour flow information are kept, such that 
the unweighted events are ready for parton shower and hadronisation. With such 
a generator at hand, a large set of studies will become available. 



4.3. Other theoretical advances (by Jean-Philippe Lansberg) 

On top of the theoretical advances mentioned above, several interesting theo- 
retical results have been obtained in recent years. Let us review some of the most 
significant ones briefiy. 

Last year, Collins and Qiu [86] showed that in general the fc-r-factorisation 
theorem does not hold in production of high-transverse-momentum particles in 
hadron-coUision processes, and therefore also for ip and T. This is unfortunate 
since many studies had been carried out successfully using fcy-factorisation (see 
references in section 3.3 of [2]), predicting mostly longitudinal yields and smaller 
CO LDMEs, in better agreement with the idea of LDME universality. 

On the side of NRQCD, Nayak, Qiu and Sterman provided an up-to-date 
proof [87] of NRQCD factorisation holding true at any order in v in the gluon- 
fragmentation channels. They showed that improved definitions of NRQCD matrix 
elements were to be used, but that this was not to affect phenomenological studies. 

Besides, the c- and 6-fragmentation approximation was shown to fail for the Prp 
ranges accessible in experiments for quarkonium hadroproduction [30] . By studying 
the entire set of diagrams contributing to tp and T production in association with 
a heavy-quark pair of the same flavour, it was shown that the full contribution 
was signiflcantly above (typically of a factor of 3) that obtained in the fragmen- 
tation approximation. The latter holds (at 10% accuracy, say) only at very large 
Pt- Pt > 60 GeV/c for ip and Pt > 100 GeV/c for T. Note that the same obser- 
vation was previously made for the process 77 — > J/ipcc [88] and also for the B* 
hadroproduction, for which it was noticed that the fragmentation approximation 
was not reliable at the Tevatron [89, 90]. 

Moreover, still in double-heavy-quark-pair production, the notion of colour- 
transfer enhancement was introduced by Nayak, Qiu and Sterman [91]. If three 
out of the four heavy quarks are produced with similar velocities, then there is 
the possibility that colour exchange within this 3-quark system could turn a CO 
conflgurations into CS ones, thus effectively increase the rate of production of CS 
pairs. They flnally discussed the introduction of speciflc new 3-quark operators of 
NRQCD necessary to deal with such an issue. 
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5. PERSPECTIVES FOR SOME (NEW) OBSERVABLES 

Standard quarkonium measurements with general purpose detectors at the LHC are 
generally related to kinematical distributions of the quarkonium decay products, 
such as differential cross section and polarisation measurements, and focus on 
decays into muons. Although these provide useful information, it is important to 
investigate the use of additional observables. In particular, it could be helpful to 
take into account not only the kinematics of the quarkonium itself, but also that 
of particles produced in association, or their nature as for instance in the study 
the production of quarkonia in association with a heavy-quark pair. 



5.1. Hadronic activity around the quarkonium 

(by A.C. Kraan) 

Here, we investigate observables that are sensitive to the hadronic activity di- 
rectly around the produced quarkonium [92]. This allows to extract information 
about the radiation emitted off the coloured heavy-quark pair during the produc- 
tion, and thereby about its production mechanism itself. To study the sensitivity of 
a typical multi-purpose LHC detector, we generated J /ip and T events in PYTHIA 
8 [93] in four production toy models: colour-singlet and three colour-octet models 
with a varying amount of shower evolution of the coloured QQ-stdXe. 

Because at lower Pr/ , shower activity is small in general, differences manifest 

themselves only at higher values of P-p , about 20 GeV/c. It must be noted that 
the energy of the surrounding particles associated with J/ ^-production is small 
(order GeV), and as such it is not a priori clear whether there is any sensitivity 
at all over the underlying-event background. Also, a careful understanding of the 
detector is necessary, so we would not recommend this analysis to be done with 
early LHC data. 

In Fig. 1 (left) we display the transverse momentum density dP^/dfi/j for J/ip^s 
between 20 and 40 GeV/c after reconstruction, in a cone around the J ftp of certain 
size R = A/(Ar7p~+~(A0)2, where 

^paround (^ j^) ^ paround ^j^ ^ ^j^/2) - P|>-ound (j^ _ ^j^/g) 

dn^ ~ 7r[(R + dR/2)2 - (R - dR/2)2] ' ^' 

Here P^'"°"°'^(R) is the sum of the transverse momentum of all charged particles 
(with Pt > 0.9 GeV/c) inside the cone of size R. In Fig. 1 (right) we display this 
variable for prompt T(lS')-events. For T-events the activity is lower than that for 
J/ip-events because the bb is heavier than the cc state, and thus has a smaller 
shower evolution. 

For prompt J/ip^s, the main background comes from non-prompt J/ip produc- 
tion. Whether one can subtract the background component from data in order to 
obtain the hadronic activity related to prompt production remains to be seen. For 
T's, the situation is easier: the background can be studied accurately by studying 
the events in the side-bands. 
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FIGURE 1. Variable dP^'''"'"'^{R)/dnu at reconstruction level for J/^P (left) and T(IS') 
(right) events. 



5.2. Associated production channels (by Jean-Philippe Lansberg) 

Another very valuable observable, likely to test the many production models 
available [1, 2], is the study of associated production channels, first in pp collisions, 
then in pA and AA. By associated production channels, we refer to ip + cc and 
T + bb. 

A first motivation for such studies is simple: similar studies carried at B- 
factories showed an amazingly large fraction of J/ip production in association with 



aje-* 



(T(e-t 



*J/i)+cc) 



another cc pair. Indeed, Belle collaboration first found [94] 

0.59^q;}5 ±0.12. Thereafter, the analysis was improved and they obtained [95] 



»j/V+x) 



to be 



a[e~^e 



J/ij + cc) 



a{e+e'^J/ij + X) 



0.82±0.15±0.14, 



(4) 



> 0.48 at 95% CL. 



Whether or not such a high fraction holds for hadroproduction as well, is a 
question which remains unanswered. Analyses at the Tevatron (CDF and D0) and 
at RHIC (PHENIX and STAR) are already possible. As computed in [30] for the 
RUN2 at the Tevatron at v^= 1.96 TeV, the integrated cross-section are significant 
(see Fig. 2 for the J/ip): 



a^J/ip + cc) X 
a{T + bb)x 



1 nb 
Ipb 



(5) 



Without taking into account any modifications of the CO LDMEs induced by 
the QCD corrections mentioned in Section 4.1, the integrated cross sections were 
found in [33] to be dominated by the CS part, similarly to the differential cross 
section in Pt up to at least 5 GeV/c for ip and 10 GeV/c for T. In other words, such 
observables can be thought of as a test of the CS contribution, for the first time 
since the idea that CO transitions would be the dominant mechanism responsible 
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FIGURE 2. Differential cross section for pp —^ Jj-tp + cc as function of the J/-0 transverse 
momentum Pt at the Tevatron {^/s= 1.96 TeV). 



for quarkonium production at high transverse momentum. If the effect of CO 
transitions is confirmed to be negigible for the T, the T produced in association 
with a bb pair are predicted to be strictly unpolarised, for any Pt (see Fig. 3). 
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FIGURE 3. Polarisation of an T produced in association with a bh pair at the Tevatron for 
^=1.96 TeV for \y\ < 0.6. 

Beside the property of discriminating between the CO and the CS transitions, 
the yield of ip in association with cc should be a priori less sensitive to the Xc 
feeddown and B feeddown (and T with bb insensitive to the Xb feeddown). Indeed, 
being suppressed by the relative velocity, the P-wave yield is expected to be smaller 
than the CS S-waves^. Note that the situation is completely at variance with the 




To be complete, let us mention the possibility to produce Xc + cc via the process gg -^ gg for 
which the two final-state gluons split into a cc pair, one of them habronising into a Xc via the 
CO mechanism. This contribution is certainly suppressed up to Pt — 20 GeV. For larger Pt, a 
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inclusive case for which it is nearly always easier to produce a P-wave than a i/j 
since this requires one less gluon attached to the heavy-quark loop. Concerning the 
ip feeddown from B, the same is expected to occur: there are only prices to pay to 
produce a ip via a B with a cc (heavier quark mass and then decay into a tp) while 
no gain in the Pt dependence since both -B(— * ^-^) + cc and tp + cc cross sections 
scale like P^^. 

Let us also mention that associated production has also been studied in direct 77 
collisions in Ultra-peripheral collision (UPC) [96]. At least for direct 77 collisions, 
associated production is the dominant contribution to the inclusive rate for Pt > 
2 GeV/c. 

To conclude, let us mention that studies can be carried on by detecting either the 
"near" or "away" heavy-quark with respect to the quarkonia. There are of course 
different way to detect the D, B, or a 6-jet, ranging from the use of a displaced 
vertex to the detection of their decay in e or /i. This has to be considered by 
also taking into account the different backgrounds. In any case, we hope that such 
measurements would provide with clear information on the mechanism at work in 
quarkonium production. 



5.3. Exclusive quarkonium photoproduction in proton and 

nucleus colliders (by David d'Enterria) 

A significant fraction of proton-proton and ion-ion collisions at collider ener- 
gies involve "ultraperipheral" electromagnetic interactions characterised (in the 
Weizsacker- William equivalent-photon-approximation [97]) by the exchange of a 
quasi-real photon. Exclusive quarkonium photoproduction at proton or nucleus 
colliders - i.e. processes of the type 'fh^Vh where V = J/ip,T and the hadron 
h (which can be a proton or a nucleus A) remains intact - has been measured 
at RHIC [98] and the Tevatron [99] and will be measured at the LHC in both 
pp [100] and pPb and PbPb [101] collisions. Exclusive QQ photoproduction offers 
an attractive opportunity to constrain the low-x gluon density at moderate virtu- 
alities, since in such processes the gluon couples directly to the c or 6 quarks (see 
Fig. 4) and the cross section is proportional to the gluon density squared (see [102] 
and refs. therein). The mass of the QQ vector meson introduces a relatively large 
scale, amenable to a perturbative QCD (pQCD) treatment. In the case of nuclei, 
the information provided by such processes is especially important since the gluon 
density is very poorly known at low-x and there are not many experimental 
handles to measure it in a "clean" environment [103]. 

The CDF collaboration has recently reported preliminary measurements of 
exclusive J/ip, ip' and T photoproduction in the dimuon decay channel in p-p 



dedicated calculation is needed. However, this mechansim would be very easily disentangled from 
the CS contributions since both c quarks are necessarily emitted back to back to the Xc and thus 
to the J/i^- 
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FIGURE 4. Schematic diagram for exclusive difFractive quarkonium photoproduction in ^A, 
7P collisions collisions. 



collisions at 1.96 TeV [99]. The cc states can be very well observed above a small 
dimuon continuum from 77 interactions (Fig. 5, left). The data is compared to 
photon-pomeron predictions as implemented in the starlight Monte Carlo [104] 
in order to try to pinpoint a possible excess which could be indicative of photon- 
odderon J/ip production. At higher masses the T(IS') and T{2S) are also clearly 
visible. Similar simulation studies in pp collisions at the LHC have been carried 
out by ALICE [105] and CMS [100]. 

In Ultra-Peripheral Collisions (UPCs) of heavy-ions the maximum photon ener- 
gies attainable are uJmax ~ 3 GeV (100 GeV) at RHIC (LHC). Correspondingly, the 
maximum photon-nucleus cm. energies are of the order WJ^^^ ~ 35 GeV (1 TeV) 
at RHIC (LHC). Thus, in 7^4 — >■ J / il) (T) A^*"* processes, the gluon distribution 
can be probed at values as low as x = My/W'^y<^ ~ 10~^(10~^). At low-x, gluon 
saturation effects are expected to reveal themselves through strong suppression 
of hard-exclusive diffraction relative to leading-twist shadowing [106]. While this 
suppression may be beyond the kinematics achievable for J /ip photoproduction 
in UPCs at RHIC, x ^ 0.01 and Ql^ ^ M^/4 ^ 3 GeV^, it could be important in 
UPCs at the LHC [101]. 

The PHENIX experiment has measured J /ip photoproduction at mid-rapidity in 
Au-Au UPCs at ^s^;^ = 200 GeV in the dielectron channel [98]. Within the (still 
large) experimental errors, the preliminary J/ip cross-section of da / dy\\y\^Q,^ = 
48±14(stat) ± 16(syst) fib is consistent with various theoretical predictions [104, 
107, 108, 109] (Fig. 5, right). The band covered by the FGS predictions includes 
the J/ijj cross sections with and without gluon shadowing [107]. Unfortunately, the 
current statistical uncertainties preclude yet any detailed conclusion regarding the 
nuclear gluon distribution. 

At the LHC energies, the cross section for T(15') photoproduction in UPC PbPb 
^t y'Sjvjv — 5-5 TeV is of the order of 150 fih [104, 110]. Inclusion of leading- 
twist shadowing effects in the nuclear PDFs reduces the yield by up to a factor 
of two, o"x = 78 /ib [110]. Even larger reductions are expected in calculations 
including gluon-saturation (Colour Glass Condensate) effects [HI]. Full simulation 
studies of input distributions generated with the starlight MC [104] have shown 
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that ALICE [105] and CMS [22] can measure well J/ip — ^ e+e^, yU+yU^ and T — > 
e+e~, /i^yU~ respectively (in different pseudorapidity ranges) in UPCs tagged with 
neutrons detected in the ZDCs. 



5.4. Bottomonium dissociation at the LHC (by David Blaschke) 

A lot of progresses have been made in the experimental and theoretical inves- 
tigations of charmonium production in heavy-ion collisions, see [112] for a recent 
review. Although CERN-SPS data from the NA50 and NA60 collaborations have 
provided a clear evidence for an anomalous J/ip suppression pattern, an unam- 
biguous theoretical approach to the phenomenon is still lacking. As elements for 
such an approach are being assembled, it becomes clear that theoretical models 
have to be supplemented by measurements of nonperturbative phenomena such as 
to provide baselines for new physics, i.e. for the studies of heavy quarkonia in hot 
dense matter as produced in high-energy AA collisions. 

While the situation with charmonium production becomes more and more com- 
plex and puzzling when comparing SPS and RHIC experiments due to intricate 
intertwining of nonperturbative initial-state, formation-stage and final-state effects 
(see Section 3.2), the bottomonium spectroscopy accessible in the upcoming LHC 
experiments may offer a much cleaner probe of the physics of dense matter. This is 
due to the following specificities of the heavier bb system relative to the cc system: 
(i) absence of down-feeding from a heavier quarkonium system, (ii) negligible regen- 
eration of bottom quarks from the thermal medium, (iii) dominance of Coulombic 
part of the bb interaction and partonic medium, so that plasma screening and ther- 
mal dissociation effects shall be under better control, (iv) both low-lying states, 
T(15') and T'{2S) are bound states in the temperature range T = 1...2 Tc, and 
are thus well-separated from the hadronisation stage. 

Here we would like to review some first baseline estimates of bottomonium 
spectroscopy for the discussion of upcoming experiments at the LHC and their 
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impact on the discussion of quark-gluon plasma (QGP) properties. Due to the 
above characteristics, bottomonium production at the LHC will allow to test the 
kinetic approach to heavy-quarkonium dissociation by thermal activation in a quark 
plasma. In this approach the survival probability is given by 



S-qi 



SrG SqgP Snuc 



exp - / TaciT)^ exp - / Tqgp{T)^ exp{-nN(J^bsL) .(6) 



where the quarkonium dissociation rate (inverse lifetime) in a QGP is determined 
by the thermally averaged breakup cross sections by (massless) quark and gluon 
impact 



QGP 



(T) 



^QGP(^) 



^ 27r2 



1=1,9 



ijj dua 



{QQ)A^Fre\ni[uj) 



(7) 



The dominant medium effect on the dissociation rate is given by the temper- 
ature dependence of the binding energy, i.e. the energy necessary to excite the 
quarkonium state to the continuum threshold where it can perform rearrangement 
reactions to open flavor states without energy cost. These inputs are provided 
from solutions of the heavy-quarkonium Schrodinger equation with a temperature- 
dependent heavy-quark potential, see Figs. 6 and 7. While the binding energies 
for J/iIj and T' drop to zero for temperatures just above Tc, the T is a bound 
state up to at least 2.2 Tc- The absolute value of the binding energy is below the 
thermal energy of the impacting partons ~ T and therefore the T production shall 
be dominated by thermal dissociation processes. 




ro 



FIGURE 6. The colour singlet QQ free energy F{r,T) vs. r at different T [f f3, f f4]. 

For the importance of the in-medium modification of the threshold, we want to 
refer to processes with quark impact, described by the Bethe-Born model (BBM) or 
the string-flip model (SFM), see [116] and references therein, results are shown in 
Fig. 8, left panel. Corresponding results for gluon-induced processes have recently 
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a function of the temperature: BBM and SFM approach give a similar cross section enhancement 
due to the lowering of the breakup threshold, from Ref. [116]. Right panel: Survival probability for 
T in a longitudinally expanding gluon plasma as a function of the plasma lifetime, from Ref. [4]. 



been discussed in [117]. The right panel of Fig. 8 gives a rough estimate at which 
level T suppression by thermal dissociation in a QGP is to be expected. A more 
detailed discussion of bottomonium dissociation at RHIC and the LHC can be 
found in [118]. 

Summarizing, we want to emphasise the role of precise bottomonium spec- 
troscopy in AA collisions at the LHC for the diagnostics of QGP properties, such 
as temperature and lifetime. A key role for the extraction of these properties do 
play microscopic theories for dissociation reactions of the states of the hh spectrum 
with their kinematic dependences. 
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6. BRIDGING THE GAP BETWEEN pp AND AA 

(by Joseph Cugnon) 

This section is devoted to a few remarks on the actual understanding of quarkonium 
production in heavy-ion colhsions in terms of the present knowledge of the prop- 
erties of quarkonium production in free space, of the properties of the quarkonium 
and on the heavy-ion collision dynamics. 

Of course, this survey starts with the status of the theoretical description of the 
quarkonium production in NN collisions. The traditional idea of a gluon fusion 
followed by the colour-octet mechanism seems to fail to reproduce the CDF data. 
However, this idea within a colour-singlet picture has been reconciled recently with 
the data for J /il) at low Pt [28] . In a simple language, the interaction between the 
c and c quarks has been added. 

As said in Section 3.3, the conventional belief is that, in heavy-ion collisions, 
quarkonium states are produced in the first A^A^ collisions, similarly as in free 
space. This is however probably not true. At high energy, the incoming heavy-ions 
should be regarded as fluxes of partons. The latter may not be involved only in 
the first collisions. The fluxes at the moment of interaction may be altered. These 
"initial state effects" on the parton distributions have been invoked to describe the 
y and Pt J /ip distribution in heavy-ion collisions. It seems that these distributions 
can be described by advocating either the cc interaction or the broadening of the 
kr distribution of the initial gluons [69]. 

Let us now examine the status of the so-called J /ip suppression. The original 
idea by Matsui and Satz [6] is that the screening of the colour forces inside a 
plasma leads to the eventual dissociation of quarkonium states. Expecting such a 
plasma in heavy-ion collisions at sufficiently large energy, it was predicted that the 
J /if) yield should be reduced. Since the original yield is not known, the idea is to 
look at the relative yield as a function of the path length in the plasma (in the 
famous Raa versus L plot). And indeed, such a suppression was observed in the 
NA38 experiment soon after. In the meantime, Hiifner and collaborators showed 
that such a suppression is obtained in a simple multiple A^A^ collision picture, 
provided an inelastic J/'i/'-nucleon inelastic cross section of a few mb is used [119]. 
This result shed some trouble for a while, but nobody believes in this scenario 
any more. Evidences for J /%p suppression has been accumulated, even if they are 
not always consistent. This is not considered as a proof of the existence of the 
plasma, for several reasons. First, the energy spectra of the produced particles 
do not show temperatures above Tc (see the interesting discussion of H. Satz in 
these proceedings [120]). The alternative probes of the plasma do not give clear-cut 
answers. Furthermore, the J/ip suppression can be accounted for thanks to several 
alternative scenarios (comovers, meson gas or fluid). 

Interesting developments have taken place in the recent years. There is more and 
more evidence that quarkonia survive in the plasma up to temperatures of the order 
of 2Tc [70], indicating that the transition would not be a pure second order phase 
transition but rather of the Kosterlitz-Thouless type [71]. Therefore, the coupling of 
a quarkonium state and the plasma has to be re-examined. The J/^/'-gluon inelastic 
cross section has been re-evaluated recently to take into account this partial 
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screening [72, 121]. Furthermore, the higher-mass resonances are not automatically 
dissolved. They should increase the J/?/' yield. Therefore it has been suggested 
that the observed J /ip suppression is largely coming from the disappearance of 
these resonances (the so-called "direct" suppression). The evidence of that is 
controversial and there is little hope that experiments at the LHC will clarify 
the issue (see Section 4.3). 

In the recent years, attention has been put on open c and open h production 
through charmed and bottomed mesons and/or jets driven by c or 6 quarks. In 
particular, the Pr-dependence indicates a strong suppression, contrasting with the 
small suppression of quarkonium (even no suppression at large Pt'-, incidentally, the 
interest in the Pr-dependence of J /if) suppression has been so luckily revived, this 
variable being less ambiguous than the L variable or equivalent). This has raised 
a strong interest in the energy loss of a heavy quark in the plasma. According to 
some authors [69, 122, 123], the radiative energy loss in colour fields may lead to 
an upper bound of the final parton energy, expressed as 

Abound /~\ JP'i T ^ ' 

where m is the mass of the quark, L is the path length and F is the colour force. 
This results from the so-called AdS/CFT correspondence and crucially depends 
upon the fact that the specific energy loss is proportional to the square of the 
energy in such an approach. From reasonable estimates, the previous equation 
leads to 

^6ow(GeV) < -f^ (9) 

where ^=1 for c and 14 for h quarks. If this is true, heavy-quark production will 
be limited to surface emission only. In addition, recent measurements indicate that 
c and h quarks participate to the fiow, which even complicates the picture of the 
propagation inside the plasma. Definitely, this issue warrants further investigation. 
In conclusion, if there has been undeniable theoretical progress on various aspects 
of the production of quarkonium in the recent years, one has to admit that our 
understanding is still far from complete. 

7. CONCLUSION 

On the verge of the start-up of the LHC, an overview of the current knowledge 
of heavy-quarkonium production reveals a situation that is somehow still not 
satisfactory. In proton-proton collisions, the interest of this process stems from the 
simple observation that the rather large scale introduced by the heavy-quark mass 
allows a separation between perturbative and non-perturbative physics, opening 
up the way to a first principle description. However, given the present data, the 
effectiveness of such an approach is yet to be confirmed. 

Concerning nucleus-nucleus collisions, the status of J/ip suppression as an indi- 
cator of the formation of the quark-gluon plasma (QGP) is also becoming more 
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shaky. There has been a great effort to use proton-nucleus colhsions to evaluate 
the so-called Cold Nuclear Matter (CNM) effects on the J /ip suppression in the 
recent years. 

Yet the situation appeared more complicated than expected concerning the nu- 
clear parton distribution functions (nPDFs) at low x and their dynamical evolution. 
Although novel experimental observables, like e.g. quarkonium photoproduction in 
ultra-peripheral nucleus-nucleus collisions, provide a new handle to shed light on 
these issues. The influence of the QGP on the suppression is not really transparent: 
the understanding of the screening of colour charges above the critical temperature, 
of the coupling of the heavy quarkonium to the plasma and of its time evolution 
seems to be less solid than a few years ago. 

On the other hand, several aspects are promising, most of them being related 
to the advent of the LHC beams. As for heavy-quarkonium production in proton- 
proton collisions, the possibility of bottomonium production, of measurements of 
polarisation and associated production of J/ip and of T with a heavy-quark pair will 
severely constrain the current theoretical models. As for proton-nucleus collisions 
the analysis of the dAu data at RHIC and the set-up of an elaborated pA program 
are certainly of importance and promising of future programs. 

Finally, for nucleus-nucleus collisions, the LHC will offer new possibilities: higher 
temperatures of the plasma and production of T states, which are expected to 
clarify the analysis. Theoretical advances concerning the interaction of the heavy 
quarkonia with the hot plasma, the slowing down of heavy quarks and the space- 
time evolution of the plasma are expected. Those advances are anyhow necessary. 
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